Shear-transformation-zone theory of plastic deformation near the glass transition 
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The shear-transformation-zone (STZ) theory of plastic deformation in glass-forming materials 
is reformulated in light of recent progress in understanding the roles played the effective disorder 
temperature and entropy flow in nonequilibrium situations. A distinction between fast and slow 
internal state variables reduces the theory to just two coupled equations of motion, one describing 
the plastic response to applied stresses, and the other the dynamics of the effective temperature. 
The analysis leading to these equations contains, as a byproduct, a fundamental reinterpretation 
of the dynamic yield stress in amorphous materials. In order to put all these concepts together 
in a realistic context, the paper concludes with a reexamination of the experimentally observed 
rheological behavior of a bulk metallic glass. That reexamination serves as a test of the STZ 
dynamics, confirming that system parameters obtained from steady-state properties such as the 
viscosity can be used to predict transient behaviors. 



I. INTRODUCTION 

For over a decade, my coworkers and I have been devel- 
oping a shear-transformation-zone (STZ) theory of plas- 
tic deformation in noncrystalline solids, [ll, H, H, 0, 0] 
Our goal has been to construct a phenomenological de- 
scription of amorphous plasticity based on physical prin- 
ciples and molecular models, and yet simple enough to 
be useful for predicting the performance of real materials. 
At the molecular level, amorphous solids are structurally 
no more complicated than fluids. They do, of course, ex- 
hibit highly non-fluidlike properties - rigidity, jamming, 
and the like. Nevertheless, their underlying simplicity 
implies that their behaviors might exhibit some degree 
of universality. My purpose in this paper is to move the 
STZ theory further toward such a description. In par- 
ticular, I propose to start with a few basic features of a 
flow-defect theory of amorphous plasticity and then to 
see how far it is possible to go using primarily symmetry, 
conservation laws, and thermodynamics. 

From its inception the STZ theory was intended 
to be an extension of the flow-defect theories of Turn- 
bull, Cohen, Spaepen, Argon and others [sl. [ol. [Tol [Til [T^ 
in which localized deformable clusters of molecules al- 
low noncrystalline solids to undergo irreversible shear 
strains in response to applied stresses. The reformula- 
tion of the STZ theory presented here is motivated pri- 
marily by the emergence of the effective disorder temper- 
ature as a key internal state variable in nonequilibrium 
theories of amorphous plasticity. The dynamic role of 
the effective temperature has been explored in remark- 
able detail recently by Haxton and Liu in their ex- 
tensive molecular-dynamics simulations of a simple, two- 
dimensional glass forming material. Much of the content 
of the present paper was developed during the attempt 
by Manning and myself to interpret the Haxton-Liu 
data. 

The relation between the population of STZ's (or flow 
defects) and an intensive variable such as the effective 
temperature has a long history. Earlier investigators, no- 
tably Cohen and TurnbuU [Ts'l and Spaepen , described 



the intrinsically disordered state of noncrystalline mate- 
rials by a free volume Vf. Those authors perceptively 
recognized that the relevant definition of w/ is not as 
an extensive excess volume measured from some densely 
packed state, but as an intensive quantity - the inverse of 
the derivative of the configurational entropy (i.e. the en- 
tropy associated with molecular conflgurations, without 
kinetic contributions) with respect to the volume. Thus 
they proposed that the density of flow defects be propor- 
tional to a Boltzmann-like factor, exp (— constant/?;^), 
and not just to Vf itself. 

In I argued that the appropriate generalization of 
free volume in plasticity theory is an effective tempera- 
ture Tg// that characterizes the state of configurational 
disorder in the system. In analogy to Vf, T^ff is the 
inverse of the derivative of the configurational entropy 
with respect to configurational energy. Tg// equilibrates 
to the ambient temperature T at high T, but may fall 
out of equilibrium at low T where disorder is generated 
by the molecular rearrangements that accompany me- 
chanical deformation. Throughout this paper, as in 
I define Teff = Tz Xi where Ez — ksTz \s s. character- 
istic STZ formation energy, so that the STZ density is 
proportional to exp(— 1/x). This is a direct analog of 
the free- volume formula and, in fact, reduces to it in the 
case of a system under constant pressure with a positive 
"effective" thermal expansion coefficient. 

With this definition of x, the "bottom line" of the anal- 
ysis to be presented here is that the STZ theory ~ in most 
but not all circumstances - can be reduced to an expres- 
sion for the plastic part of the rate-of-deformation tensor. 
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(1.1) 



supplemented by an equation of motion for x, which 
appears here only in the Boltzmann-like prefactor. In 
Eq. (jl.ip , s is the deviatoric stress tensor and T is the bath 
temperature. The tensor function fij{s,T) contains all 
the dynamical details pertaining to STZ transitions; the 
prefactor exp(— 1/x) determines the density of STZ's. 
This clean separation between the s and x dependence 
of the plastic strain rate was emphasized by Shi et al. 
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[T^ as a characteristic feature of the STZ theory. It is 
derived here in Section IIIII and the equation of motion 
for X is derived in Section [VII 

This paper starts, in Sec[TTl with a reexamination of 
the basic assumptions upon which the STZ theory is con- 
structed. The STZ equations of motion are derived in 
Section Hill and a general expression for the yield stress 
is derived and interpreted in Section [TVl The latter re- 
sult brings new physical insight to the choice of material- 
specific quantities such as the STZ transition rates. Sec- 
tions |V] and IVH respectively, contain a derivation of the 
STZ formula for the Newtonian viscosity and a discussion 
of effective-temperature dynamics. As a way of putting 
all these concepts together in a realistic context, the main 
part of the paper concludes in Sec. IVIII with a reexam- 
ination of the extensive set of rheolo gica l measurements 
of a bulk metallic glass by Lu et. al., [17|]. That reexami- 
nation, a reworking of ^ and ^Tsll , serves as a test of the 
STZ dynamics, confirming that system parameters ob- 
tained from steady-state properties such as the viscosity 
can be used to predict transient behaviors. Finally, Sec- 
tion IVIlIl contains a few brief remarks about outstanding 
problems. 



II. BASIC ASSUMPTIONS OF THE STZ 
THEORY 

The starting point for all of the analysis in this paper is 
a dynamic model of a noncrystalline material in which a 
disordered arrangement of molecules (or bubbles, or sand 
grains, etc.) interact with each other via short ranged 
forces. This system is subject to driving forces that cause 
it to deform and, in some circumstances, to flow contin- 
uously. It is dynamic in the sense that it possesses an 
intrinsic time scale, say tq, that characterizes the rate 
at which it responds to microscopic perturbations. In 
molecular materials, tq is an internal vibration period of 
order femtoseconds. The theoretical challenge is to un- 
derstand the mechanism by which these fast, molecular 
interactions generate viscous responses on time scales of 
order seconds or longer. 

To construct a theory of amorphous plasticity for such 
a model, I have found it useful to watch computer- 
generated, moving pictures of its behavior as it under- 
goes slow, steady, shear deformation. These pictures can 
be found on the websites of the authors of Refs. [l^ and 
[T9| . Their most striking feature is that the irreversible 
molecular rerrangements are sporadic, short lived, spa- 
tially isolated events. These are the STZ transitions. 
They occur within a persistently noisy environment in 
which thermally and mechanically generated fluctuations 
bring groups of molecules into and out of conflgurations 
in which they can undergo irreversible shear transfor- 
mations. That is, the fluctuations create and annihi- 
late STZ's. At low temperatures, below the glass transi- 
tion, rearrangements are driven primarily by the applied 
stress; and an STZ, once formed, rapidly undergoes a 



shear transition if it is aligned favorably with respect 
to the stress. At higher temperatures and small driving 
forces, i.e. in the viscous regime above the glass transi- 
tion, STZ transitions are thermally driven and the ap- 
plied stress simply biases them the direction of the aver- 
age plastic shear rate. In either limit, the STZ transitions 
are rare inelastic events occurring within an otherwise 
solidlike elastic material. 

It is especially easy to observe STZ rearrangements in 
two-dimensional models, where their cores are Tl events 
in which two nearest-neighbor molecules move away from 
each other and two next- nearest neighbors come together. 
No two STZ transitions ever seem to look exactly the 
same. The Tl events at their cores occur at varying 
angles with respect to the shear direction; and varying 
numbers of nearby molecules participate in the overall 
motions. This statistical variability of the STZ's is clearly 
due to the fluctuating environment in which they are be- 
ing created and annihilated. 

Experimental observations of memory effects in sim- 
ple amorphous materials tell us that the STZ's cannot 
be structureless objects. The Bauschinger effect is one 
example where the system remembers the direction in 
which it has been deformed, and responds differently ~ 
more compliantly or less so - to further loading in differ- 
ent directions. The natural way to include such effects 
in the theory is to let the STZ's possess internal degrees 
of freedom that carry information from one event to the 
next. The simplest such possibility is to assume that, 
during their lifetimes, the STZ's are dynamic, anisotropic 
objects whose populations and orientations are deter- 
mined by the loading history. 

Once a shear transformation occurs in some direction, 
the molecules that composed the STZ resist further shear 
in the original direction but may be especially suscepti- 
ble to a reverse shear. That is, the flrst transition re- 
distributes the local stresses in such a way as to favor a 
reverse transition if the stress changes sign. There is no 
strong requirement that the reverse transition bring the 
molecules back to exactly their original positions; but it is 
this approximate picture that suggests a two-state model 
of STZ's. The two-state behavior has been observed di- 
rectly by Lundberg et al. [l^ in parallel bubble-raft ex- 
periments and two-dimensional foam simulations. These 
authors find that a Tl event is often mechanically re- 
versible (although still dissipative) if the direction of the 
applied shear is reversed shortly after the event occurs. 
On the other hand, if the shear is not reversed until after 
other events have occurred nearby, then the memory is 
lost, i.e. the STZ is annihilated. 

Thus, the two-state STZ's are the memory carriers 
in amorphous materials. If they are randomly oriented 
and present in substantial numbers, the material is de- 
formable in all directions because there are STZ's avail- 
able to respond to all directions of applied stress. If they 
are rare, deformation is difficult in any direction. If they 
are all oriented in the direction of a shear stress, then 
further deformation in that direction is impossible - the 
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material is jammed - but response to a reverse stress 
is relatively easy. Annihilation of existing STZ's and 
creation of new ones without orientational bias are the 
mechanisms by which memory is lost in these systems. 

During its lifetime, each STZ has two preferred spa- 
tial orientations. One of these orientations, say the "+" 
state, will be closest to the direction of the deviatoric 
stress tensor s and therefore most stable. In two dimen- 
sions, the other STZ state, say "— ", is perpendicular 
to the first and is therefore least stable with respect to 
s but most stable with respect to — s. The situation is 
more complicated in three dimensions where the two STZ 
states need not be symmetrically oriented with respect to 
s; but it must still be true that, except in very special 
cases, if we are given a stress tensor s and a pair of STZ 
orientations, we can identify the "±" orientations unam- 
biguously. 

Let the symbol w denote the orientation of an STZ, 
that is, the orientation of the axes along which its ± 
states are defined, and let the symbol i denote other 
characteristics of the STZ such as its size or its transition 
threshold. Then let n±{u},i) be the associated STZ den- 
sity. In accord with the preceding discussion, a master 
equation for this set of these densities has the form 

Toh±{uj,£) =TZi{±s^)n^{uj,£) - Tlii^s^) n±{uj,£) 



-[r(s)+p(T)] 



(2.1) 



where tq is the molecular time scale mentioned previ- 
ously. The right-hand side of Ea. (|2.ip consists of two 
parts, one that couples the STZ's to the stress in a way 
that conserves the total number of STZ's and the config- 
urational entropy, and a second that describes the rate 
at which STZ's are created and annihilated and therefore 
governs the entropy flow in the system. 

The first two, entropy conserving terms on the right- 
hand side of Eq. (|2.ip describe the rates at which the 
STZ's transform between their and "— " states. 

T^e{suj) is the rate of transitions from " to and 

is the projection of s onto the axes defined by u, . By 
symmetry, the rate of reverse transitions is 'R-i{—s^). 

The second pair of terms on the right-hand side de- 
scribes the rate of STZ creation and annihilation. Here, 
(r -|- p)/tq is an attempt frequency, i.e. a noise strength, 
consisting of incoherent mechanical and thermal parts, 
r(s) and p(T) respectively. The exponential function 
exp (—Oi/x) is the a priori probability of occurence of 
a state with energy Eg — ksTz Of, and rioo is a reference 
density of order an inverse molecular volume. r(s) is a 
measure of the strength of the mechanically generated 
noise that accompanies plastic deformation; it vanishes 
when the driving force s is absent. An explicit expression 
for r(s) will be derived in Sec. IIIII 

p{T) is the super-Arrhenius part of the rate of ther- 
mally activated molecular rearrangements. Its molecular 
origin lies at the heart of any theory of the glass tran- 
sition, and is discussed in detail - albeit quite specula- 
tively " in my paper on an "excitation-chain theory of 



glass dynamics. [20, l2l| p{T) vanishes when T is less 
than the glass transition temperature Tq and is equal to 
unity when T is above the super-Arrhenius region, say, 
T > Ta- It is convenient to write 



P{T) 



(T) 





for T > To 
for T < To 



(2.2) 



where a{T) vanishes for T > Ta and diverges, perhaps 
like (T — To)^^, as T approaches Tq from above. When 
there is no external driving and the effective temperature 
Teff = xT^z equilibrates to T, the STZ creation rate is 
proportional to exp [—Ee/kB T — a{T)], which is the "a" 
relaxation rate. Importantly, the activation energy ap- 
pearing here is the sum of Arrhenius and super-Arrhenius 
parts. Eg and kBTa{T) respectively. 

Two aspects of Ea. (|2.ip require extra attention. First, 
note that r(s) and p{T) are assumed to be independent 
of the STZ label £; they are noise strengths that ap- 
ply equally to all molecular rearrangements. Second, the 
creation and annihilation part of Eq. (j2.1D is written as a 
single detailed-balance relation in which the equilibrium 
distribution always is proportional to the Boltzmann fac- 
tor exp {—Oi/x), with x rather than T in the exponent. 
The a priori probability of forming a configurational de- 
fect always is determined by the effective disorder tem- 
perature. 



III. STZ EQUATIONS OF MOTION 
A. State Variables and Time Scales 

The next order of business is to extract as much phys- 
ical insight as possible from the structure of the STZ 
theory summarized by Ea. (|2.ip . making approximations 
that retain this structure, but without being specific 
about the model-dependent ingredients of the rate fac- 
tor 7^. 

Because only the STZ's couple to the applied stress, 
the plastic rate of deformation tensor Z?fj can be written 
in the form 



X (^d^jito) [Uiis^) n^{Lo,£) - Uei-s^) n+(co,£)] ^ .(3.1) 

Here, eo is a shear increment of order unity. The angular 
brackets denote an average over STZ orientations u con- 
sistent with the "±" constraints (see [Hj for details) . The 
symbol / d£ denotes a weighted sum over the other STZ 
properties including transition thresholds. The traceless, 
symmetric tensor dij (w) projects these transitions onto 
the i, j axes. 

For two dimensional systems, Pechenik [5] showed that 



(3.2) 



where e(w) is a unit vector at an angle uj relative to 
a principal axis of the stress, with — 7r/4 < uj < 7r/4. 
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If this principal stress axis is at an angle, say (p, with 
respect to the x axis, then s,y — sdij{4>), where s is 
the (signed) deviatoric stress. Throughout the following 
discussion, I adopt the convention that positive values of 
s denote stresses in the direction that drives "— " to "+" 
transitions. 

The situation is more complicated in three dimensions 
where the tensor dij{uj) must contain more than just di- 
rectional information. Note, for example, that we need 
two stress values plus three angles to specify the devia- 
toric stress tensor s. Since the tensorial versions of the 
STZ equations to be used here must look the same in 
three as in two dimensions, I omit the detailed three di- 
mensional analysis. 

In principle, we should solve Eq. (|2.1|) separately for the 
n± (w, at each uj and each I, and then use those n± (cj, H.) 
to evaluate Z?fj in Eq. (|3.ip . The results of such a calcu- 
lation may be important in some situations. For exam- 
ple, the distribution over orientations lo may be quantita- 
tively relevant when the system is driven strongly away 
from equilibrium or when the orientation of the stress 
changes abruptly in time. Similarly, at small stresses, 
the sum over STZ thresholds implied by / M det ermines 
the extent of plastic deformation; only the STZ's with 
low thresholds undergo transitions before the system be- 
comes jammed. Such calculations, however, would be 
more laborious than is necessary for most practical pur- 
poses. Note especially that only a narrow range of STZ 
thresholds is likely to be dynamically relevant. STZ's 
with anomalously low thresholds will have large forma- 
tion energies and therefore be very rare; that is, they will 
be suppressed by the weight factor implicit in J M. On 
the other hand, STZ's with high transition thresholds will 
not contribute appreciably to the deformation rate; they 
will be suppressed by the rate factor Tli{s^). Thus the 
sum over £ seems likely to be dominated by STZ's with 
a single characteristic formation energy that already has 
been denoted by ksTz- 

If only a narrow range of values of £ is dynamically rel- 
evant, then the integration over € in Eq. (|3.ip contributes 
just a numerical factor that we can assume already has 
been incorporated into the reference density Uoo- Then, 
dropping the variables £, rewrite Eq. (|3.ip in the form 



The total STZ density, normalized to noo, is 



(3.3) 



(3.4) 



which case we can remove this sum from inside the inte- 
gration. As a result, the first term on the right-hand side 
of Eq. (|3.3p can be written in the form 



2 no. 



= eoA(d 



,(cj) (7^(s^)-7^(-s^) 



(3.5) 



Similarly, we may assume that the sum TKsS) -f TZ{—Si^) 
depends only weakly on the orientation and can be re- 
placed by an w-independent function of the magnitude of 
the stress, specifically. 



7^(s^)+7^(-s^) 



7^(s)-|-7^(-s) =C{s). (3.6) 



where s — (1/2) Sij and TZ{±s) denotes the value 
of TZ for a stress of magnitude s oriented along the "±" 
directions. Then the second term on the right-hand side 
of Eq. (|3.3p can be written in the form 



2n~ 



(w) [(7^(s^) + 7^(-s^)) {n+{uj) - n_(co))] ^ 



eo AC(s) m,y. 



(3.7) 



where 



n-|_(ct;) — n-(Lu) 
n+iuj) + n-{uj) 



(3.8) 



is a deviatoric tensor that describes the average STZ ori- 
entation. In analogy to Eq. (|3.7p . rewrite Eq. (|3.5p in the 
form 

eoA(d„(c^) (7^(s^)-7^(-s„))) -eoAC(s)^T(s), 



(3.9) 



where 



^T(.-) 
s 



s 



n{s) - TZj-s) 

^^(iyr^^Pi) 



(3.10) 



With these definitions, the plastic rate of deformation 
tensor in Eq. p.ip becomes 



ToDf^^eoC(s)A 



fn-s) 



rrii 



(3.11) 



Then return to the master equation for the STZ densities, 
Eq. p.ip . to deduce equations of motion for nnj and A. 
Our assumption about a narrow range of relevant values 
of £, e.g. thresholds, means that we write 0^ = 1, so 
that all the relevant STZ formation energies are of order 
Ez = Tz as anticipated in the definition of x- Using 
the preceding definitions of rriij and A, we find 



which defines the dimensionless density A. The sum 
n+iuj) + n-{uj), as opposed to either of these terms sep- 
arately, ought to be approximately independent of w, in 



2C{s) -^T(s) 



!!!y_!_e-i/x. 



A 



(3.12) 
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and 



where 



ToA = r* 



(e-i/'^-A); (3.13) 



r*°* = r(s) + p(r). 



(3.14) 



The variables A and rriij , and the roles that they play 
in Eqs. (l3.11l - I3.13|) . immediately tell us a great deal 
about the nature of the STZ theory. A is the fraction 
of molecular sites occupied by STZ's. To be consistent 
with the assumption of dilute, weakly interacting STZ's, 
A must be small, no larger than 10"'^, and usually very 
much smaller. The deviatoric tensor rriij is the average 
STZ orientation. By definition, its magnitude is less than 



Pechenik argued that the symmetry preserving and 
physically intuitive way to determine T is to assume that 
it is proportional to the rate at which the work of defor- 
mation is dissipated irreversibly, and thus is converted 
into the disordered configurational fluctuations that cre- 
ate and annihilate STZ's. This assumption determines T 
uniquely if one invokes the second law of thermodynamics 
by requiring that this dissipation rate be non-negative. 
The argument starts by writing the first law of thermo- 
dynamics in the form (using summation convention): 



D 



-C(s)A ^T(s) 

To 

d 



— *(A,my) + Q(sy,A,my). (3.15) 
The left-hand side of Eq. (|3.15p is the rate per unit volume 



or equal to unity. In Eq. (|3.1ip . m.ij plays the role of a at which plastic work is done by the stress Sij. On the 

right-hand side, "if is the recoverable, state-dependent, 
internal energy density associated with the STZ's. ^ 
must be proportional to the density of STZ's; therefore 
it is convenient to write it in the form: 



back stress, consistent with the idea that STZ's already 
aligned in the direction of the stress impede further de- 
formation in that direction. The fact that the rate factor 
TZ{s) is a non-negative, monotonically increasing function 
of its argument means that T{s) is a monotonic function 
of s that vanishes when s — 0. Its magnitude, like that 
of rriij ) is bounded by unity. 

We also deduce from Eas. (|3.11l - I3.13P that there are 
two qualitatively different time scales in the STZ theory. 
The plastic strain rate determined by Eq. (|3.1ip is a rate 
per unit volume; it is proportional to the small quantity 
A because it scales with the density of STZ's. No pref- 
actors A appear in the expressions for rhij and A on the 
right-hand sides of Eqs. p.l2p and (|3.13p . These equa- 
tions describe how individual STZ's respond to changes 
in their environments, and the fact that A is missing as 
a prefactor in those equations implies that rriij and A re- 
spond to perturbations much more rapidly than does the 
rate of plastic deformation. (The factor exp (— l/x)/A on 
the right-hand side of Eq. (|3.12p is of order unity accord- 
ing to Eq. (|3.13p .) We have not yet written an equation of 
motion for the effective temperature x; but it should be 
clear that x is a slow variable. (See Sec. IVIH The time 
derivative x must be proportional to the rate per unit vol- 
ume at which configurational entropy is generated during 
plastic deformation; thus, like the expression for x 
must contain a prefactor A. 



B. The rate factor T 

Determining the rate factor F that first appears here in 
Ea. (|2.ip has been one of the most challenging problems 
in the STZ theory. This factor has been chosen incor- 
rectly in much of the earlier literature in the field. In 
our original STZ paper |lj], Falk and I used an expres- 
sion for F that turns out to be correct at sufficiently low 
temperatures; but we did not propose a systematic ratio- 
nale for it. The problem was solved by Pechenik [3, Q, 
who showed that F can be determined by thermodynamic 
arguments alone. 



*(A,my) = eoAV'(m). 



(3.16) 



Because ■0 is a scalar function of just , it must depend 
only on TO = (1/2) mijuiij. Also assume that the ef- 
fective temperature x is so slowly varying compared to A 
and rriij that it may for the moment be taken to be a con- 
stant, and need not be included as an explicit argument 
of the internal-energy function ^E". 

The last term on the right-hand side of Ea. (|3.15p . 
i.e. Q, is the energy dissipation rate per unit volume. 
Pechenik's hypothesis is that F is proportional to the 
rate of energy dissipation per STZ. That is, 

Q(sy, A, TOy) = So — AF(sy , A, rriij). (3.17) 
To 

where sq is an as yet undetermined factor with the di- 
mensions of stress. 

The next step is to use Eqs. (l3.12p and (|3.13p to eval- 
uate the time derivatives in Eq. (|3.15p and solve the re- 
sulting equation for F or, more conveniently, F*°*. The 
result is: 



F*°*(s,j,A,my) = 



^ s 



1 



A(A,to) 



— [fn) 

So m 



-sop{T) 



(3.18) 



where 



-i/x 



-i/x 



A(A, to) = So — TO 'ip\m) — h Tp{ra) ^ — -j^ 1 

(3.19) 

The quantities F, F*°* and p(T) all must be non-negative 
- the first because of the second law of thermodynamics, 
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the second and third because they are non-negative rate 
factors. The condition that the numerator in Eq. (|3.18p 
remains non-negative for all values of is sufficient to 
determine ip(m) up to an (unnecessary) additive con- 
stant. To find ip{rh), compute the inverse function of 
T; specifically, find the function 



So 



(3.20) 



(sjj/so) — {'^ij I "rn) S,{ffi) change sign at the same value 
of Sij , and therefore that the product of these factors is 
never negative. 

At this stage in the development, it is useful to take 
advantage of the fact that A relaxes so rapidly that we 
may replace it in the preceding formulas by its steady- 
state value, A exp(— 1/%), and let 



A(A, fh) ^ sq {1 - m £,{fn)). 



(3.21) 



By definition, T(s) vanishes at s = and goes smoothly 
to 1 as s ^ oo. Therefore fh^(rh) vanishes like fh^ at 
m = and diverges at to = 1. So long as to remains 
in the range, < m < m^aa; , whe re fhmax ^{mmax) = 
1, the denominator A in Eq. (j3.18p remains positive as 
required. Moreover, the dynamics of the system never 
allows an initially small to to reach ffimax because the 
dissipation rate diverges at that point. 

With these simplifications, the equation of motion for 
TOij, Eq. (|3.12p . becomes 



1 



To rriij 
C 



A(A, m) 
(,-) (^r(s-) - TO, 



So 



myp(r)) .(3.22 



As mentioned previously, this is a stiff differential equa- 
tion; there is no factor A — exp(— 1/x) on the right-hand 
side to produce slow relaxation. Therefore, for most pur- 
poses, we can set rriij = and replace to^ elsewhere by 
its stress dependent value for which the right-hand side 
of Eq. (|3.22p vanishes. Moreover, for an isotropic system 
in which the only orientations are set by the stress tensor 
s;, , we must have 



to. 



■M(s) 



(3.23) 



where M{s) is the stationary solution of Eq. (|3.22p with 
the preceding ansatz. Specifically, 



M{s) 
2s' 



so_ 
2s 



i + ^r(s-)-,^' 

So 2C(s) 



(l + -r(s-) 
So 



P{T) 
2C(s) 



-4 — r(s). (3.24) 

So 



We then have 
roi?^'(s-,x) = eoC(s)e-i/^ '-f [T{s) - M{s)\, (3.25) 

which has the form anticipated in Eq. (|l.ip . 
For later reference: 

r.o.(.) _ 2 C [r - Af] [(s7so) - e(M)] + P{T) ^^^^^ 



and 



such that T[so^(to)] = to. Then, because T is a mono- 
tonically increasing function of its argument, the choice 
Tp'^fh) = ^(to) ensures that both {sij /s)T{s) — rriij and r(s) — 



1 - Mi{M) 



2C[T- M] [(s/so) - ^ (M)) + M e(M) p{T) 



l-Af^(M) 

(3.27) 

where C, T, and M are all understood to be functions 
only of s and T . The last expression reduces to the for- 
mula for r(s) postulated in py, but it does so only for 
temperatures low enough that p{T) = (i.e. below the 
glass transition), and that reverse STZ transitions are 
neghgible, 7^(-s) « 0. Then r(s) ^ 1 and £,{M) « 0, so 
that 



So^e-i/>^r(s-)«7^^'s,,. 



To 



(3.28) 



This relation between the STZ production rate and the 
applied power density has been confirmed recently by 
Heggen et al. [1^ in the context of conventional flow- 
defect theories. 



IV. THE DYNAMIC YIELD STRESS 

One of the most notable features of the STZ theory 
is the natural way in which a yield stress emerges from 
its basic structure. To see this, consider the case where 
the bath temperature T is below the glass transition, 
so that the spontaneous STZ annihilation and creation 
rates, proportional to p(T), both vanish. This condition 
does not necessarily mean that the transition rate TZ{—s) 
vanishes, as assumed in the strictly athermal version of 
the STZ theory, [i, Thermal fluctuations may induce 
backward STZ transitions even when the system is in a 
completely inviscid glassy state, i.e. when it has infi- 
nite linear shear viscosity - which is the only situation in 
which there is a sharply deflned yield stress. 

Setting p{T) = in Eq. p.22p . or simply evaluating 
Eq. ([XM)l in this limit, we flnd M{s) to be 



M{s) 



T{s) for s < Sy, 
sq/s for s > Sy, 



(4.1) 



where Sy is the value of s at which the two branches of 
M{s) cross; that is, Sy is the solution of 

SyT{sy) = so. (4.2) 

Clearly, Sy is a dynamic yield stress. When s < Sy, 
Eqs.dSSSl) and (|iT|) imply that Dfj = 0; the system is 
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jammed and no plastic deformation is occurring. When 
s > Sy, D^j is nonzero. The two cases of the function 
M(s) shown in Ea. (|4.1[) are the stable branches of the 
steady-state solutions of Eq. (|3.22p ; thus the yield stress 
occurs at an exchange of stability. Note that the de- 
nominator on the right-hand side of Ea. (|3.22p vanishes 
only at TO = So /sy = T{sy), that is, at the value of to 
corresponding to the exchange of stability at s = Sy] 
therefore to is dynamically constrained to remain be- 
tween these limits as anticipated in the discussion pre- 
ceding Eq. (13221) • 

As defined here, the yield stress Sy is an intrinsic, 
steady-state property of a material; it does not depend on 
the history of deformation or even on the material's ini- 
tial state of disorder. This quantity can appropriately be 
called the "ultimate yield stress" or the "minimum flow 
stress." It is quite different from the "peak stress," which 
will appear in the transient stress-strain curves computed 
in Sec lVIII and which has no such intrinsic meaning. 

Equation (|4.2p unambiguously determines the propor- 
tionality factor So that appears in the Pechenik relation, 
Ea. p.l7p . in terms of the yield stress Sy. Because Sy 
generally depends on temperature, pressure, and per- 
haps other state variables, so must also depend on those 
variables. The temperature dependence of so may be 
especially important near a glass transition, where ther- 
mal fluctuations become increasingly effective in assisting 
molecular motions over energy thresholds, so that Sy is a 
decreasing function of increasing T. In fact, sq is a more 
fundamental quantity than Sy. It remains well defined 
above the glass transition, where Eq. (l4.2p is no longer 
valid and where, strictly speaking, there is no yield stress, 
but a function such as £'fj(s, x) in Eq. (|3.25p may still ex- 
hibit visible stress dependence as s increases through sp- 
in the literature on flow-defect theories of amorphous 
plasticity, it is generally assumed that - in the language 
of the STZ theory - the yield stress is determined by 
some intrinsic feature of the transition rate TZ{s). Sup- 
posedly, 7?.(s) rises sharply when the stress grows large 
enough to drive molecules in an STZ over an energy bar- 
rier, and the barrier stress is the yield stress Sy. This 
assumption cannot be correct. For example, in the low- 
temperature limit mentioned at the end of Sec lIIIi where 
7^(-s) « , Eg. iprU)) tells us that T{s) « 1, so that 
Sy ~ So, independent of whatever function we choose for 
TZ{s) and without any necessary relation to the dynamic 
ingredients of the rate factor. 

How, then, are we to understand so? The easy answer 
is that we obtain it from experiment or calculate it in 
some approximate way, and then build it into our choice 
of TZ{s). The chosen Tl{s) determines the behavior of 
^ij i^' x) I and thus fixes the relation between stress and 
plastic flow at all stresses, including near so- But, had 
we chosen a "yield stress" different from so in our choice 
of TZ{s), we would not have changed Eq. (|4.2p very much 
or obtained a significantly different value of Sy. In fact, 
we do not know whether any characteristic stress appear- 
ing in TZ{s) is related to so- Therefore, simply inserting 



a physically motivated choice of sq into both 7l{s) and 
Eq. (|4.2p . while perhaps sensible, evades the basic ques- 
tion. In principle, we need independent estimates of both 
So and TZ{s). 

To find an independent estimate of sq, return to the 
original Pechenik relation, Eq. (|3.17p . The left-hand side 
of this equation is Q = T dS/dt, where dS/dt is the rate 
per unit volume at which the entropy S is increasing as 
a result of plastic deformation. On the right-hand side, 
writing A = cxp(— identify 



fo 
To 



-i/xr(s) 



dSc 
~dt 



(4.3) 



where [dSd dt)„iech is the rate at which configurational 
entropy in the form of STZ's is being produced during 



mechanical deformation, and VL 



is the volume per 



molecule. The factor kBVz is an approximation for the 
entropy of an STZ, which means that vz is roughly the 
number of molecules in an STZ. Thus Eq. (|3.17p takes the 
form 



dS 



dt kBVz 



So 51 f dSc 
'dt 



(4.4) 



From this we deduce that the fraction of the total dis- 
sipation rate Q that produces STZ-like configurational 
disorder is vz ks T/sq Q. In other words, an amorphous 
material yields more easily if it converts a larger fraction 
of the work of deformation into configurational disorder 
instead of ordinary heat. 

The relation between sq and entropy generation in 
Eq. (|4.4p is not purely formal; it has physical content, 
and therefore should be useful for evaluating sp. To il- 
lustrate this possibility, we can make rough estimates of 
both sides of Eq. (|4.4p . For this purpose, assume that the 
temperature is very small, well below the glass transi- 
tion. On the left-hand side of Eq. (14. 4p , consider any ir- 
reversible molecular rearrangement driven by an applied 
shear stress. The system first deforms elastically as the 
molecules are driven to an unstable threshold, and then 
all of this elastic energy is converted into some combina- 
tion of kinetic energy (heat) and configurational energy 
as the system relaxes to its rearranged state of equilib- 
rium. The local shear strain at threshold must be of 
order unity; thus we have TdS/dt nj, where is the 
shear modulus, and the total shear rate 7 is a measure 
of the rate at which these rearrangements are occurring. 
To evaluate the right-hand side of Eq. (|4.4p . note that 
this term is equal to the rate at which STZ's are be- 
ing created, each adding an entropy increment of order 
ks per molecule. Again, assume that the relevant low- 
temperature rate of events is proportional to 7, so that 

n {dSc/ dt)mech ~ ks Vz 7- 

With these estimates, Eq. (l4.4p becomes simply 



So 



Vz' 



(4.5) 



This is the relationship found by Johnson and Samwer 
[23| . who pointed out that the low-temperature yield 
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strain Sy//i is an almost universal constant, of order 0.03, 
for thirty different bulk metallic glasses. A literal inter- 
pretation of Eq. (|4.5p would imply that vz ~ 30, which 
may be reasonable for these complex, multicomponent 
materials. However, Eq. (|4.5p is not much more than a 
dimensional analysis, and should not be taken so seri- 
ously. A complete calculation would include a theory of 
the relation between Q and the elastic driving forces, and 
would require a better estimate of the relation between 
the size of an STZ and its entropy. Both of those calcu- 
lations are well beyond the scope of this paper. 



rearrangements and not by the STZ transition rate it- 
self. In that connection, note also that I am mak- 
ing a different assumption here than the one I used 
in 4], where I guessed that the rate factor 7?.(s) van- 
ished super-Arrheniusly at the glass transition, and that 
p{T)/C{Q) ^ 1 at that point. That assumption was in- 
consistent with the fact that, even though their linear 
viscosities vanish, glassy materials do undergo plastic de- 
formation at low temperatures and large driving forces. 
(See the discussion at the beginning of Section IVlin 



V. NEWTONIAN LINEAR VISCOSITY 

It is convenient at this point to use the STZ equations 
derived so far to obtain an expression for the Newto- 
nian linear viscosity. Consider the case of temperatures 
just above the glass transition, where p{T) is small but 
nonzero, and compute the plastic strain rate in the limit 
of vanishingly small driving stress s. The quantity M{s) 
defined in Ea. (|3.24p must vanish linearly as s — > 0, so 
the term ifis/so in the numerator of Eq. (|3.22p may be 
neglected and, with Eq. (|3.23p . the equation for M(s) be- 
comes 



Therefore, 



2C 



M{s) 



)-M]=Mp{T). 



T{-s) 



1 



P(T) 
2C(s) 



P{T) 
2C(0) 



1 



(5.1) 



(5.2) 



For steady-state motion at temperatures above the glass 
transition, in the limit of infinitesimally small deforma- 
tion rate, the effective temperature T^ff becomes the bath 
temperature T, so that x T/Tz- Using Eqs. ((3?25)) . 
(|5.ip . and (|5.2p . we find (for pure shear) 



DP'- 



-DP' 
yy 



2to 1 



PiT) 
2C(0) 



The Newtonian viscosity is 



?7jv(r) = lim 



onP' 



To 

eor'(O) 



2C(0) 



To 



eor'(O) p(T) 



= ?7o exp 



^ + «(r) 



(5.3) 
(5.4) 



where ry^^ = (eo/To) '?''(0), and the super- Arrhenius 
bmction Q!(T) is defined in Eq. (|2.2p . The approximation 
in Eq. (|5.4p is valid only at temperatures close enough to 
the glass transition that p(T) ^ C(0), and thus neglects 
a potentially important temperature dependence of the 
prefactor 779 at larger T . 

Note that, especially in its limiting form near the 
glass transition, the viscosity is determined almost en- 
tirely by the super- Arrhenius rate of thermally activated 



VI. EFFECTIVE-TEMPERATURE DYNAMICS 

We now need an equation of motion for the dimen- 
sionless effective temperature x- There is an extensive 
literature on this subject. References that I have found 
useful include [H H [II, [13, [11 . Much of what follows is 
based directly on work by Liu and colleagues, especially 

As in [4'] , the equation of motion for x to be used here 
is an approximate statement of heat balance. It has the 
form 



CeffTzX - Teff 



dSc 
~dt 

dSc 
~dt 



mech 



the 



1 



X 



xTz 



T 



(6.1) 



On the left-hand side, Ce// is an effective specific heat 
per unit volume of the form Ceff = where co is 

a dimensionless number of order unity. Thus this term is 
a rough estimate of the rate at which the configurational 
heat content is changing as a function of time. 

The two terms on the right-hand side of Eq. (|6.ip are 
the rates at which the configurational heat content is be- 
ing changed, respectively, by mechanical work and by 
thermal fluctuations. The factor (dSc/dt)mech, intro- 
duced in Eq. (|4.3p , is the rate at which conflgurational en- 
tropy is produced in the form of STZ's during mechanical 
deformation. The definition of the effective temperature 
implies that we multiply this rate by T^ff to obtain the 
heat production. In the second term, {dSc/dt)therm de- 
notes the rate of configurational entropy production by 
thermally induced molecular rearrangements. This term 
is closely related to the thermal STZ creation rate, pro- 
portional to the function p(T) introduced in Eq. (|2.ip : but 
it requires further consideration. 

Both the factors (1 — xlx) and (1 ~ xTz/T) appear- 
ing in Eq. (j6.ip introduce physical mechanisms that are 
beyond the basic STZ assumptions. The second is the 
easiest to understand. It says simply that, at tempera- 
tures above the glass transition, the effective temperature 
X relaxes toward T/Tz] and it models this aging effect 
by a conventional linear law of cooling. The first fac- 
tor, i.e. the modification of the mechanically driven rate 
of entropy production, is the more interesting. It says 
that, at low temperatures where thermal fiuctuations are 
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negligible, x relaxes to a steady-state value denoted here 
by x(q), and it assumes that this relaxation - like the 
cooling law used in the thermal term - is linear. 

To interpret the mechanical factor, think about bath 
temperatures below the glass transition, so that molec- 
ular rearrangements are not spontaneously activated by 
thermal fluctuations but must be driven by externally 
applied forces. Those forces, in effect, "stir" the system 
at a rate, say, 7, which may be the norm of the rate-of- 
deformation tensor. The only intrinsic time scale in the 
system at low temperatures is tq; therefore, define the 
dimcnsionless stirring rate to be g = 7x9. Consider first 
the limit q <^1. If, in order to achieve steady-state sta- 
tistical equilibrium, each molecule must have changed its 
neighbors at least once or twice, then the equilibration 
time is irrelevant; only the magnitude of the deformation 
makes a difference. After long, slow stirring, the state of 
the system is characterized by an effective disorder tem- 
perature that is independent of the stirring rate or any 
details of the stirring mechanism. Denote this value of x 
by Yojthat is, x ^ Xo as 9 ^ 0. (In earlier papers, e.g. 
BSlll, Xo was denoted by Xoo-) 

More generally, if q is not negligibly small, then the 
steady-state value of x is x{l)- The g-dependence of x 
is likely to be nontrivial for foams or granular materials, 
where tq is not microscopically small, and experimentally 
accessible strain rates may be comparable to Tq~^. This 
situation also may arise in amorphous molecular mate- 
rials when strain rates are very large, for example, near 
crack tips or at the centers of shear bands. As in p^). 
the q dependence of x may be quite interesting. 

The linear approximations made in both of the terms 
on the right-hand side of Eq. (|6.ip imply that this equa- 
tion cannot be used to describe very large excursions from 
steady-state equilibria, x cannot be far from x{q) when 
{dSc/dt)mech is large and {d Sc/dt) therm is small; nor can 
X be very far from T/Tz in the opposite situation. 

Return now to the entropy-production terms in 
Eq. (pT|) . The factor {d Sc/dt)raech as given in Eq. (|i3)) 
appears here, because the mechanical generation of con- 
figurational entropy requires the STZ mechanism. The 
problem of evaluating the thermal factor {d Sc/dt)therm 
is more problematic and interesting. By definition, the 
effective temperature pertains to all the configurational 
degrees of freedom of the system, not just those that cou- 
ple to external shear stresses. Above the glass transition, 
thermal fluctuations may generate a variety of different 
kinds of defects, with formation energies different from 
ksTz, and these will contribute to the configurational 
entropy. With this possibility in mind, but keeping as 
close an anology to Eq. (|4.3p as possible, I propose that 
the thermal factor have the form: 



dSc 
~dt 



therm 



To il 



P{T) 



(6.2) 



Here, I have replaced the mechanical factor T{s) in 
Ea. (|4.3p by the thermal factor p{T) on the assumption 
that these functions play the same roles as noise strengths 



as they did in Ea. (|2.ip . To account for the possibility 
that non-STZ configurational disorder may be included 
in the entropy production, I have replaced exp (— 1/x) by 
K exp (— /3/x); where P Tz is an activation free energy 
for defect formation, and k is a dimcnsionless rescaling 
of the factor vz/Vl. Note that the parameters fi and k 
play the same roles here as they did in [J; but now it 
is more obvious why (3 should be smaller than unity, be- 
cause it is a rough approximation for a free energy that 
may include more disordered states than those described 
just by STZ's. In Q, I showed that systems with /3 < 1 
may be especially susceptible to shear-band formation. 
With the preceding assumptions, Eq. (|6.ip becomes 



tqCoX 



-i/xr(s)x 
e-f'/^p{T)x 



1 



X 



xiq) 

TzX 



T 



(6.3) 



where cq = (cq/coI'z)- r(s) is shown in Eq. (|3.27p . This 
form of the x equation differs slightly from its original 
version, Eq.(3.5) in Importantly, use of Eq. (l4.3p sim- 
plifies the final result and eliminates the need for an extra 
approximation in the derivation. Note that x{q) natu- 
rally appears in the denominator inside the first brack- 
eted term on the right-hand side of Eq. (|6.3p . Thus, large 
values of x do not produce unphysical behavior as would 
have happened previously. 

The effective temperature is emerging as a remarkably 
powerful concept for understanding the nonequilibrium 
properties of amorphous materials. It plays a key role in 
the analysis of large-scale deformation of metallic glasses, 
first published in [4| , and reworked here in Section rVIII In 
that case, the dynamics of x as determined by Eq. (|6.3p 
control the plastic response of the system to transient 
changes in external driving. I know of no other mech- 
anism capable of quantitatively explaining the observed 
relaxation phenomena. 

Some of the most interesting recent developments in- 
dicate that, in steady-state nonequilibrium situations, x 
can be used "quasithermodynamically" as if it were an 
ordinary temperature. That is, x can be used as an inde- 
pendent intensive variable, along with the thermal tem- 
perature, the pressure, etc., in equations of state for ex- 
tensive quantities such as the volume or internal energy. 

Perhaps the first example of a quasithcrmodynamic 
analysis is that of Shi et al. [l6j , who postulated a linear 
equation of state for the potential energy as a function 
of X and used this to interpret their molecular-dynamics 
simulations of shear banding. Manning, Carlson, and I 
[S^ l used Eq. (|6.3p and the data in [16] to show that shear 
banding occurs here via a nonlinear, transient instability. 
Another quasithcrmodynamic analysis appears in a pa- 
per by Bouchbinder, Procaccia and myself Q, in which 
we used an athermal version of STZ theory to interpret 
molecular dynamics simulations of amorphous silicon by 
Demkowicz and Argon. [10, Hi], HI, [s^ The latter authors 
measured not just stress-strain curves but also the frac- 
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tion of atoms that were in liquidlike, as opposed to solid- 
like, nearest-neighbor configurations. We found that this 
liquidlike fraction, in steady-state nonequilibrium situa- 
tions, obeys a quasithermodynamic equation of state as 
a function of x- With this understanding, we were able 
to account quantitatively for the simulation results, in- 
cluding the time dependent transients in the liquidlike 
fraction that Demkowicz and Argon observed near the 
onset of loading. 

The most strikingly unexpected and speculative qu- 
asithermodynamic role played by x is suggested by 
the molecular dynamics simulations of Haxton and 
Liu.[l^ These authors simulated a simply sheared, two- 
dimensional, glass-forming material over three decades of 
steady-state strain rates 7 and for bath temperatures T 
ranging from about one tenth of the glass transition tem- 
perature To to about twice Tq. By measuring pressure 
fluctuations, they independently determined values of x 
for each value of 7 and T. For T < Tq, Manning and I 14 1 
propose that their observed relation between 7 and x is 
a direct analog of the relation between the a relaxation 
rate and the bath temperature T near a conventional 
glass transition. More specifically, the relation between 
the dimensionless strain rate q and the low-temperature 
(p = 0), steady-state x = xil) in Eq. (|6.3p can be written 
in the form 



1 1 

— exp 
90 



-T + aeffix) 



(6.4) 



which is a direct analogy to the final form of the Newto- 
nian viscosity in Eq. (|5.4p . Here, q^^ and A are constants 
analogous to ijo and Tz- ae//(x) is a super- Arrhenius 
function that has the same form as a{T) in Eqs. (|2.2[ ) 
and ()5.4p . diverging at an effective Kauzmann tempera- 
ture xo and vanishing above some xa- In other words. 
Manning and I interpret g to be a dimensionless rate 
at which molecular rearrangements are driven by x fluc- 
tuations below Tq. According to this extreme version 
of the quasithermodynamic hypothesis, q should exhibit 
the same behavior - including a glass transition at xo - 
as that which occurs for thermally driven rearrangements 
above Tq. Another interesting prediction of this hypothe- 
sis is that X diverges at a finite value of the dimensionless 
strain rate, i.e. at g = go- If true, the material would 
"melt," and the solidlike STZ assumptions would fail at 
a well defined upper limit of the driving strength. 



VII. STZ THEORY OF DEFORMATION IN A 
BULK METALLIC GLASS 

To date, the most complete applications of the STZ 
theory in the interpretation of laboratory experiments 
have been our analyses 1, [11] of the metallic glass data 
published by Lu et al. [13| . These papers were written be- 
fore we understood several crucial aspects of the theory. 
For example, we had not recognized that the separation 
of time scales discussed in Section IIIII implies that the 



parameter eo never appears by itself in the equations of 
motion, but always in the combination cq/tq. As a result, 
the values of eo and tq cited in those papers were dramat- 
ically incorrect, although the errors mostly cancelled out 
in the final results. In the present interpretation, eo ~ Ij 
and To is of order femtoseconds; thus cq/to ^ lO^^sec"^. 

Another example is that, in those earlier papers, we 
had not understood the relationship - or lack thereof 
- between the STZ transition rate 7?.(s) and the yield 
stress, and therefore were not thinking carefully enough 
about the physical basis of TZ{s). We can now write an 
expression for this rate factor that is directly related to 
molecular mechanisms, without being constrained to in- 
clude the yield stress explicitly in that formula. More- 
over, the excitation-chain theory of anomalously slow re- 
laxation near the glass transition [l^, |2l[ implies that the 
super- Arrhenius factor p{T) determines the rates of spon- 
taneous creation and annihilation of STZ's and thermal- 
ization of the effective temperature, but does not belong 
in TZ{s) itself - as was assumed incorrectly in [4j. 

It is useful, therefore, to conclude this paper by revis- 
iting the metallic glass data, both to update the earlier 
papers and to show how the theoretical ideas developed 
in this one can be brought to bear on real- world phenom- 
ena. 

Lu et al. measured deformations of the bulk metal- 
lic glass Vitreloy 1, Zr4i.2 Tiia.s CU12.5 Niio 6622.5- To 
a good approximation, their system consisted of a uni- 
form bar with uniaxial compressive stress applied, say, 
in the x direction, and stress-free surfaces normal to the 
y and z directions. This total stress tensor aij has only 
one nonzero element, axx = o", which is the experimen- 
tally reported stress. Therefore u = {^i/2)sxx — \fis. 
Similarly, the measured plastic strain rate is Dg'^, = 
(2/V3) where D^'' is the magnitude of the plastic 
rate-of-deformation tensor ZPfj in Ea. p.25p . The Newto- 
nian viscosity is 



1. ^XX 1. ^ 

rjN — um — = lim 



^0 nP' 

^^xx 



(7.1) 



therefore Eq. ()5.4p remains unchanged. 

We next must specify the STZ transition rate TZ{s). 
For the moment, return to a notation in which s = s and 
fa — m have signs, because backward transitions with 
s < play a role in the analysis, as is clear in Eqs. (|3.10p . 
In earlier versions of the STZ theory, we always chose 
the simplest possible forms of this transition rate on the 
assumption that whatever data we had available from 
experiments or simulations would not justify additional 
theoretical complications. Manning and I departed from 
this purely phenomenological app roach in |14l | , primarily 
because the Haxton-Liu data [13] that we were interpret- 
ing extended over an exceptionally wide range of stresses 
and strain rates, and therefore required us to construct 
and test a model that was more physically motivated 
than the earlier ones. 

The transition rate that Manning and I used in p^ . 
which I will adopt here, includes an Eyring-like activa- 
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tion factor at small stresses, and a smooth transition from 
Eyring to power-law behavior at large stresses. Specifi- 
cally, 



R{s) — exp 

We then have: 
C{s) = exp 

cosh 

and 



n/2 



(7.2) 



' T 



cosh(s//i) 



^ sinh(s//i) 



1/2 



T{s) = tanh 



T 



sinh(s//i) 



so that 



^ (m) — fi arcsinh 



T 



arctanh(TO) 



(7.3) 



(7.4) 



(7.5) 



The first factor on the right-hand side of Eq. (|7.2p is 
the Eyring rate in a form similar to that used in [l[, 
where the exponential function of s//2 causes the rate to 
saturate at large s. The energy ks Te is the height of 
the Eyring activation barrier; and we expect Te/Tz < 1. 
The parameter jl is the stiffness of the barrier at its peak. 
It appears here ostensibly in the same place that /I ap- 
peared in m, but now it has no direct connection to the 
yield stress. The Eyring factor explicitly expresses the 
fact that, even at temperatures below the glass transi- 
tion, the barrier opposing STZ transitions must be ap- 
preciably smaller than the STZ formation energy. In the 
limit of small stress, s jl: 



71(3) « exp 



Te 
' T 



(1 - s/fl) 



therefore 



and 



C{s) w e""^^/^ cosh 



TeS 



T (s) w tanh 



Te s 
7^ 



(7.6) 



(7.7) 



(7.8) 



The second factor on the right-hand side of Eq. ()7.2p 
converts the saturated Eyring function at large s to a 
power law: 



Cis) 



T{s) « 1. 



(7.9) 



Here, Si is a stress scale that should be of order Sj,; that 
is, the crossover between Eyring behavior and large-scale 



plasticity should occur near the yield stress. The power 
law has a physical interpretation. In dissipative systems, 
and when strain rates are not too large (7T0 <ti 1), we 
expect n = 1, corresponding to linear friction or viscos- 
ity. In the opposite situation, as in (isl [l3|. where rates 
are large, and the motion is controlled by hard-core colli- 
sions rather than the details of molecular interactions, we 
find Bagnold scaling [35] with n = 1/2. That is, if there 
is no natural energy scale in the problem, then dimen- 
sional analysis requires that the stress be proportional 
to the square of a rate; and if TZ{s)/tq is the only avail- 
able quantity with dimensions of inverse time, we have 
TZ{s) ~ s^/^. In the experiments of Lu et al, stresses 
never exceed the yield stress, and strain rates are always 
very much less than T(^^; therefore, from here on, I use 
the linear law with n = 1. 
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FIG. 1: (Color online) Viscosity of Vitreloy 1. The red cir- 
cles are the experimental measurements reported in [l3| ■ The 
solid curve is the theoretical fit described in the text following 
Eg. (001) . 

With this choice of the STZ transition rate, the New- 
tonian viscosity in Eq. (|5.4p becomes 



miT) 



ToflT 



1 



1 ^TE/T-a{T) 



For simplicity, I use the approximation 



„Tz/T+q(T) 



a(T) 



T~To 



e T^-To, 



(7.10) 



(7.11) 



which makes a smooth transition from Vogel-Fulcher be- 
havior near Tq to Arrhenius behavior above Ta- Fig[T] 
contains the viscosity data for Vitreloy 1 shown in Fig. 10 
of [l3|. As pointed out by Masuhr et al this data 
does not fit neatly into the usual Vogel-Fulcher scheme. 
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In fact, it is more easily fit by the Cohen-Grest formula 
[3^ . which contains no low temperature divergence, and 
which is what I used in However, the data can be fit 
by Eq. (|7.10p . which is consistent with other parts of the 
present analysis. 

The solid curve in Fig[T] has been plotted using 
ro/i/eo = lO^'^ Pa-sec, Tq = 250 K, Tz = 16,000/^, 
Te = 3, 000 if, Ta = 1, 000 iiT, Ti = 31,OOOi4: and 
a = 3. Note that the glass transition temperature Tq 
is very small, consistent with the Cohen-Grest analysis, 
and that Ti needs to be remarkably large in order to fit 
the data. (Note also, however, that the effective value 
of Ti in the range of interest is very much smaller be- 
cause of the exponential cutoff.) The ratio Te/Tz ^ 0.2 
seems reasonable, i.e. the Eyring barrier is substantially 
smaller than the STZ formation energy but has roughly 
the same energy scale. Similarly, if co/tq — lO^^sec"^, 
then fl ~ 10^ Pa, which is about one tenth the yield 
stress, and means that the individual STZ's are substan- 
tially but not excessively more deformable than the sys- 
tem as a whole. The effective activation energy in the 
high temperature region, where a{T) = 0, corresponds 
to Tz + Te = 19, 000 K, which is roughly consistent with 
earlier estimates. [3611 
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FIG. 2: (Color online) Tensile stress as a function of the scaled 
strain rate 2 rjN 7. The data points, with temperatures as 
indicated, are taken from Lu et al. The three solid red 

curves, from bottom to top, are theoretical predictions for 
temperatures T — 573 K, 643 K and 683 K respectively. 

Lu et al. [l3| summarize their steady state measure- 
ments of flow stress versus strain rate at eight differ- 
ent temperatures by showing that their data nearly col- 
lapse, with appreciable scatter, to a single curve of stress 
as a function of normalized strain rate 2 rj^ 7, where 
7 = DP^ = (2/V3) Z)P'. Their results are shown in FigH 
along with the STZ predictions for the lowest and high- 



est temperatures used in the experiments, T = 573 K 
and 683 K respectively, and an intermediate tempera- 
ture, T — 643 K. For computing the theoretical curves 
in Figl21 I have used Eq. (|3.25p to compute D^'' and the 
steady-state version of Ea. (|6.3p with /3 — 1 to compute 
X- The experimental values of 7 are all less than or equal 
to O.lsec"^, therefore <? = 7T0 < 10~^^ is completely 
negligible. Accordingly, in computing I have used 
x{q) « x(0) = Xo- Then, in addition to using the param- 
eters determined by fitting the viscosity 777V as described 
above, I fit the data in Fig[2]by choosing sq — 1.1 GPa 
(the measured room-temperature tensile yield stress di- 
vided by VS), si = 1.0 GPa, k = 0.5, and xo = 0.6. 

Having fixed the viscosity parameters, and having cho- 
sen not to use a temperature-dependent yield stress or to 
vary /?, I found strikingly little leeway in determining 
the remaining parameters n and xo- I did this by fitting 
the theory to the experimental points for T = 643 K, 
as shown by the middle curve and the filled triangles in 
Figl2l It then becomes clear that some of the apparent 
scatter in the experimental data is a systematic trend 
predicted by the theory. The data collapse is trivial at 
small strain rates where we are seeing just linear viscosity. 
At larger strain rates, the stress as a function of strain 
rate starts to flatten out as it approaches the yield stress, 
with - interestingly - the stresses at lower temperatures 
mostly falling below those at higher temperatures. The 
data go out only to about half the yield stress; therefore 
the crossover from small-stress to large-stress behavior, 
which is sensitive to the details of Tl{s), is not explored 
by the experiments. 

A more sensitive test of the theory is provided by the 
transient response seen by Lu et al in their constant 
strain-rate experiments. To interpret these experiments, 
we must include the elastic part of the stress-strain rela- 
tion. Assuming that the elastic and plastic parts of the 
rate-of-deformation tensor are simply additive contribu- 
tions, we can write the xx component of the equation of 
motion for the stress in the form 



1 da 



7-^5P'(s-,x), 



(7.12) 



where E is Young's modulus and a = V3s is the tensile 
stress. Then replace time t by 7 = 7i, and use Eq. (|3.25p 
to evaluate Z)p'. Equation (|7.12p becomes 



ds 

dj 



E 



2eo 
7^0 



'i/^C(s) (r(s)-M(s) 



, (7.13) 



where, according to data provided in Lu et al., E = 
i?/V3 = {E/uy) So = 50 Sq. This equation for 5(7) must 
be solved along with Eq. l|6.3p for xil), which becomes 



dx 
c?7 



1 



Co 7 To 
Kp{T) (1 



Xo 



TzX 
T 



(7.14) 
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At this point, the only free parameters are the dimen- 
sionless effective specific heat Cq and the initial values of 
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FIG. 3: (Color online) Tensile stress versus strain at three 
different temperatures as shown. The strain rate in all three 
cases is 7 = O.lsec."^ The initial effective temperatures, read- 
ing from the top curve to the bottom, are 730 K, 715 K, and 
675 K respectively. This figure is to be compared to Fig.l in 

m 

A selection of theoretical stress-strain curves, i.e. ten- 
sile stresses a as functions of strain 7, is shown in Figsl3] 
and 3) The six separate stress-strain curves shown here, 
with different temperatures and strain rates, constitute 
all of the curves shown by Lu et al in their Figs.l and 
2 for which their system is not driven so hard that it 
fails before deforming plastically. The curves shown here 
exhibit characteristic stress peaks at low temperatures 
and large strain rates. Both the locations of the peaks 
at 7 = 0.05 and their relaxation to constant flow stresses 
at about 7 = 0.2 are quantitatively consistent with the 
experimental data. The relaxation rate is slightly sen- 
sitive to the value of cq, which must be of order unity. 
The choice cq = 0.5, as used here, seems optimal. At 
low temperatures and large strain rates, the theoretical 
curves exhibit sharp cusps at their peaks, at about the 
same places where the experimental curves break off, in- 
dicating that the sample has failed. That behavior is 
illustrated in both figures by the curve for T = 643 
7 = 0.1, which is shown here as continuing to where it 
drops to a flow stress of about 0.5 GPa, although the ex- 
perimental curve is not shown in |17j as continuing much 
beyond the peak. 

The one adjustment that I have made in computing 
these stress-strain curves is in the choice of the initial 
values of the effective temperatures Xi — xil — If 
the experimental samples were equilibrated at their de- 
formation temperatures, then Xi would be equal to T /Tz] 
but that estimate produces theoretical stress peaks that 



are higher and sharper than the experimental ones, espe- 
cially at the lower temperatures and higher strain rates. 
I can correct these discrepancies without changing the 
peak positions or the rate at which the system relaxes to 
the flow stress if I adjust the x^'s to fit the peaks. My 
estimates of the Xi's, shown in the form of effective tem- 
peratures Ti = Xi^^z, are shown in the figure captions. 
Note that these initial effective temperatures remain well 
below xoTz = 960X. 
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FIG. 4: (Color online) Tensile stress versus strain at three 
different strain rates as shown. The temperature in all three 
cases is T = 643 K The initial effective temperatures Ti, read- 
ing from the top curve to the bottom, are 730 K, 725 K, 687 K 
and 650 K respectively. This figure is to be compared to Fig. 2 
in [13] 

These experimental samples were quenched to temper- 
atures below the glass temperature, where they had val- 
ues of X equal to the Active temperature at which their 
configurational degrees of freedom fell out of equilibrium 
with the heat bath during the quench. That fictive tem- 
perature seems likely to have been quite high, in most 
cases larger than the deformation temperature T, be- 
cause these samples could not have been quenched slowly 
enough to remain in thermal equilibrium very far into 
the glassy region. When the samples were then held at 
temperatures T during rheological measurements, x may 
have moved toward T, but would not have reached that 
temperature except in cases where T was large or the 
strain rate was very small. Therefore, I have adjusted 
the initial effective temperatures as shown. 

It seems to me that the success of the STZ theory 
in predicting transient behavior is a strong indication 
that the theory is capturing the underlying physics of 
dynamic plasticity. As noted above, essentially all of the 
system parameters were determined by steady-state mea- 
surements, so that there were no numbers that could be 
adjusted by orders of magnitude in the transient calcu- 
lations. That this worked accurately - despite the fact 
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that the theory relates molecular processes that occur on 
femtosecond time scales to macroscopically slow plastic 
deformations - increases my confidence in the basic fea- 
tures of the theory. Specifically, this analysis seems to 
be a stringent test of the assumption that the rates of 
change of STZ populations are determined by the rate at 
which disorder is generated during deformation, and the 
accompanying idea that disorder in amorphous materials 
is accurately described by an effective temperature. 

VIII. CONCLUDING REMARKS 

I close by stating just a few brief opinions about open 
questions and further directions for research. 

As implied in the text, I think that the least well de- 
veloped element of the theory presented in this paper is 
the derivation of the thermal entropy-production term, 
{d Sc/dt)therm, givcu in Eq. (|6.2p . This term seems likely 
to be strongly model dependent because it involves more 
than just the STZ degrees of freedom. It also may con- 
trol some especially interesting physical properties such 
as shear-band formation. 

Perhaps the most fundamental theoretical challenge is 
to understand the limits of validity of the STZ theory. 
How and when does it break down at high deformation 
rates? Or during large excursions from steady state be- 
havior where the separation of time scales described in 
Sec. mil becomes invalid? We know that, at some point, 
strongly driven amorphous materials must change from 
deforming slowly like solids to flowing rapidly like liquids. 
Where and how does this happen? Manning and I may 
have seen a clue in ^14] , where it appears that the effective 
temperature diverges at a large but finite strain rate of 
order Tq^. This line of investigation has been opened by 
the work of Haxton and Liu Il3l . More simulations and 



experiments along these lines should be very interesting. 

Looking from a broader point of view, I have argued 
elsewhere [S^ that a principal goal of research in solid 
mechanics ought to be to bridge the gap between atom- 
istic physics and engineering practice. My prime exam- 
ple of how far behind we are in this area is fracture me- 
chanics. It is well known that advancing cracks undergo 
instabilities that are qualitatively similar to sidebranch- 
ing instabilities in dendritic crystal growth. The latter 
instabilities have been well understood for about half a 
century, and that insight has provided the basis for major 
advances in solidification processing. As yet, we have no 
comparable understanding of the analogous instabilities 
in fracture. In [40| . I suggested a way in which plasticity 
theory might be brought to bear on this problem; but, at 
that time, the STZ theory was not well enough developed 
for much progress to be made. We are now beginning, 
in fs^l for example, to understand in an STZ context 
how driven amorphous materials may become unstable 
against spatially nonuniform deformations such as shear 
bands. Reference [Zl^l is explicitly an attempt to move 
in the direction of ordinary fracture. The STZ theory 
ought to be relevant at least to slow, ductile failure in 
either shearing or tensile modes. Can it also predict the 
behavior of fast brittle cracks? 
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